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Abstract 

We use the factorization method to find the exact eigenvalues and eigenfunctions for a particle 
in a box with the delta function potential V{x) = X5(x — xq). We show that the presence of the 
potential results in the discontinuity of the corresponding ladder operators. The presence of the 
delta function potential allows us to obtain the full spectrum in the first step of the factorization 
procedure even in the weak coupling limit A — > 0. 
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I. INTRODUCTION 



The time-independent form of Schrodinger's equation in the presence of a potential V 
can be written as^ 



where ip is the wave function, m is the mass of the particle, and E is the eigenenergy. 
Equation can be solved exactly only for a few potentials. A particle in a box and a 
particle in a delta function potential are two well-known and instructive examples.- The 
former can be used to describe quantum dots and quantum wells at low temperatures,- 1 ^ 
and the latter can be used as a model for atoms and molecules.- 

The solution for a particle in a box with a delta function potential has been investigated 
using a perturbative expansion in the strength of the delta function potential A.- Exact 
solutions have been obtained for the weak (A — > 0) and the strong (1/A — > 0) coupling 
limit s.- 

In this paper we discuss the solution for a particle in a box with a delta function potential 
using the factorization method and show that the presence of the delta function simplifies 
the factorization procedure. In this way we find the full spectrum of the Hamiltonian in 
the first step of the factorization method. We also show that this result applies in the weak 
coupling limit A — > 0. Note that if we put A = from the beginning, we need to continue 
the factorization procedure to find each eigenvalue in each step. 

II. PARTICLE IN A BOX WITH A DELTA FUNCTION POTENTIAL 

Consider a particle in a one-dimensional box of size a with the delta function potential 
V(x) = A S(x — Xq) = A 5(x — pa), where < p < 1. In this case Eq. ([!]) takes the form 



where ip n (x) and E n are the corresponding eigenf unctions and eigenvalues, respectively. 
Because of the boundary conditions, ?p n (x) = for x < or x > a, the form of the 
eigenfunctions inside the box is 



V 2 ^ + Vi) = Eijj, 



2m 



(1) 







(0 < x < pa) 



(3) 
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FIG. 1: First three stationary states i/j n ( x ) f° r 2m\/h 2 = 8, a = 3, A = 1, and p = 1/2. 



where fc„ = ^2mE n /h. The continuity condition of the wave function at x = pa gives 
A/B = sin[(p — l)/c n a]/sin(pfc n a). Because the delta function is infinite at x — pa, the 
first derivative of the wave function is not continuous and the relation between the left and 
right derivatives of the wave function can be obtained by integrating Eq. (J2j) over the small 
interval (xq — e, x + e) 

dip n (x), 2m f pa+e 



dx 



\pa+e 



dx 



\pa—e 



V(x)ijj n (x) dx 

' pa—e 

^^?Pn(pa). 



(4a) 
(4b) 



We substitute the eigenf unctions in Eq. ([3]) in Eq. ( 14b j) and obtain the quantization 
condition^ 



k n sm(k n a) 



2mA 
IF 



sm(pk n a) sin[(p — l)k n a}. 



(5) 



The solutions to Eq. (jSJ) give the energy spectrum of the Hamiltonian, E n = h 2 k^/2m. 
Figure [TJ shows the first three stationary states for 2m\/h 2 = 8, a = 3, A = 1, and p = 1/2. 



III. THE FACTORIZATION METHOD 



To calculate the eigenvalues and eigenf unctions of a Hamiltonian operator H, we can 
use a general operational procedure called the factorization method. In this method the 
Hamiltonian is written as the product of two ladder operators plus a constant, H = a^a + E. 
These operators are used to obtain the eigenf unctions of the Hamiltonian. In contrast to 
the simple harmonic oscillator, one ladder operator is usually not sufficient to form all the 
Hamiltonian's eigenfunctions and another ladder operator is needed for each eigenfunction. 
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This method was first introduced by Schr6dinger-~- and Dirao^ and was further developed 
by Infeld and HulUi and Green.— The spirit of the factorization method is to write the 
second-order differential operator H as the product of two first-order differential operators 
a and a', plus a real constant E. The form of these operators depends on the form of the 
potential V(x) and the factorization energy. 

The procedure for finding the ladder operators and the eigenfunctions is as follows.— We 
find operators a±, a 2 , 03, . . . and real constants Ei, E 2 , E3, . . . from the recursive relations^ 

a\ai + Ei = H, (6a) 

a\a 2 + E 2 — aia\ + Ei, (6b) 

a\a 3 + E 3 = a 2 a\ + E 2 , . . . . (6c) 

More generally 

a ] n+l a n+1 + E n+1 = a n a\ + E n , (n = 1, 2, . . .). (7) 

Also assume that there exists a null eigenfunction (root function) |£„) with zero eigenvalue 
for each a n , namely, 

a„IU = 0. (8) 

Hence, E n is the nth eigenvalue of the Hamiltonian with the corresponding eigenfunction^ 
(up to a normalization coefficient) 

\E n ) = o|4---4-ilU- (9) 

Figure [2] shows a schematic diagram describing the relation between root functions, the 
stationary states, and the ladder operators. As can be seen from Eq. (Q, a\a\ is equal to H 
except for a constant. Thus, a\ must have a linear momentum term to be consistent with 
the kinetic energy part of the Hamiltonian. According to Eq. ([7j), each of the annihilation 
operators a n should also have a a linear momentum term. Thus, a n can be written as 

[P + ifn(x)}, (10) 



y2m 

where P is the momentum operator and f n (x) is a real function of x. Although these 
operators are not hermitian (a) n = (P — if n (x)) 7^ a n ), their product is hermitian 

t 1 d2 1 1 1 2 , ^ df n . . 
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FIG. 2: A schematic diagram describing the relation between the root functions £ n , the stationary 
states ^> n , and the ladder operators a n , a^. 

We are now ready to find the ladder operators and eigenenergies for our problem. First, 
we consider Eq. ([6|) for n = 1 

a\a x + Ei = H. (12) 
From the form of the Hamiltonian ((2]) and Eq. ( fTTT) we can rewrite Eq. ( |T2l as 

2m 2m ^ 2m dx ~^ 1 2m ( ( ) 

or equivalent ly 

^-fi + £-^ + E 1 = X6(x-pa). (14) 
2m 2m ax 



Note that in contrast to Eq. (J2J), Eq. (jTjj) is a nonlinear first-order differential equation. 

To solve Eq. (fl4|) . we consider the left-hand and right-hand sides of the delta function 
potential separately. For these regions, Eq. (JH]) reduces to 

2m 2m ax 

which has the solution 

h(x) = ^E lC ot[^^^{x - &)], (16) 
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where b is a constant of integration. We require that fi(x) be finite in the range < x < a, 
where it is the solution for our problem. Because the singularities of the cotangent function 
are tt radian apart, we choose the points x = and the singularity points of f\{x). 

At these points the potential and hence the Hamiltonian are infinite and f\ (x) can be infinite 
at the boundaries. Hence, we have 



fi(x) 



(x < pa) 



j y/2m£h 

|V2r^cot[^f^(>-&)], (x>pa), 



(17) 



where (y/2mEi/fi)(a — b) — ir. To fix the value of E±, we need to use the discontinuity 
relation of the ladder operators, which can be obtained by integrating Eq. (1141) over a small 
interval (pa — e,pa + e) 

2m A 

f 1 (pa + e)-f 1 (pa-e) = —. (18) 



Equation (fl8|) shows that the presence of the delta function results in the discontinuity of 
fi (x) at x = pa. We use Eq. ffl7|) and b 



^2mE i < cot 
Because cotfvr 



■(pa 



irhy/2mEi, to rewrite Eq. fTTBl as 

2mA 



cot 



y/2mE 



a) 



h y y/2mEi 
cot a, Eq. ( JT9l) can be written in the form 

2mA 

ki {cot[(p — l)kia\ — cot(pkia)} 



-(pa) 



,2 ' 



or 



k\ sm(kia) 



2mA 



sin(pkia) sin[(p — l)k\a], 



(19) 



(20) 



(21) 



where k\ = y/2mEi/h. Equation ( |20l) is similar to Eq. ()5]). Note that, although we calculated 
the ground state energy of the Hamiltonian in the first step of the factorization method, we 
found the full spectrum. This conclusion applies to the weak coupling limit A — > 0. In this 
limit sin(/c 1 a) = from Eq. fl2TT) . which gives the full spectrum of a particle in a box with 
k\ = nir/a. This result has an interesting consequence. If we use the factorization method 
to calculate the energy levels of a particle in a box with A = 0, we cannot find all of them 
in the first step of the procedure. In fact, we should use the recursive relations to find each 
eigenvalue in each step.— Thus, contrary to the Schrodinger equation, the presence of the 
delta function in the factorization method simplifies the problem and repeated application 
of the recursion relations is not necessary. We need only to replace the index 1 with n in 
Eqs. (TTTj) and (|2Tj). 
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To obtain the root functions, we rewrite Eq. as 

— - — h ihk n cot(k n x) J £„,(x) = 0, (x < pa), (22a) 
i dx J 

+ ihk n cot[k n (x - 6)1 J £ n (x) = 0, (x > pa), (22b) 
j as / 

where £ n (x) = (x|£„). It can be easily checked that the following is a solution to Eq. f l2"2"j) 

{A sin(/c n x), (x < pa), 

(23) 
— B sin[/c n (x — 6)], (x > pa). 

Because in the first step of the factorization procedure the root function £i(x) is equal 
to the first eigenfunction ^i(x) (see Eq. ([9]) and Fig. |2]) and we have replaced the index 
1 by n, all the root functions are equal to the eigenfunctions of the Hamiltonian, that is, 
£,n{x) = ipn( x )- This result is similar to Eq. fl3]) for x < pa. Moreover, because k n b = k n a — sx, 
we have sin[/c„(x — 6)] = sin[/c n (x — a) + ir] = — sin[/c n (x — a)]. So, these eigenfunctions are 
equivalent to the solutions in Eq. ([3]). 

IV. SUGGESTED PROBLEMS 

Problem 1. Use the recursive relations ((6]) to show that \E n ) in Eq. ([9]) is the eigenfunction 
of the Hamiltonian with E n as its eigenvalue. 

Problem 2. Show that the eigenvalues E\, E2, -E3, . . . form a monotonic increasing se- 
quence, that is, Ei < E2 < E3 < 

Problem 3. Apply the factorization method to a particle in a box with no delta function 

potential. In the first step show that k\ = ix/a. Continue the factorization procedure to 

find all the eigenenergies and eigenfunctions. Hint: consider a n = [P + ic n cot(d n x)}, 

V2m 

where c n and d n are real constants. 
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